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ABSTRACT 

The high-density star formation typical of the merger /starburst events that 
power the large IR luminosities of Ultra Luminous Infrared Galaxies (ULIRGs) 
(LiR,(8-1000/im)>10 12 L ) throughout the Universe results to extraordinarily 
high cosmic ray (CR) energy densities of UcR~fewx (10 3 -10 4 ) UcR,Gai permeat- 
ing their interstellar medium (ISM), a direct consequence of the large supernovae 
remnants (SNRs) number densities in such systems. Unlike far-UV photons em- 
anating from their numerous star forming sites, these large CR energy densities 
in ULIRGs will volumetrically heat and raise the ionization fraction of dense 
(n>10 4 cm -3 ) UV-shielded gas cores throughout their compact star-forming vol- 
umes. Such conditions can turn most of the large molecular gas masses found 
in such systems and their high redshift counterparts (~10 9 -10 10 M Q ) into giant 
CR-dominated Regions (CRDRs) rather than ensembles of Photon-dominated 
Regions (PDRs) which dominate in less IR-luminous systems where star forma- 
tion and molecular gas distributions are much more extended. The molecular 
gas in CRDRs will have a minimum temperature of Tki n ~(80-160) K, and very 
high ionization fractions of x(e)>10 -6 throughout its UV-shielded dense cores, 
which in turn will fundamentally alter the initial conditions for star formation 
in such systems. Observational tests of CRDRs can be provided by high-J CO 
and 13 CO lines or multi-J transitions of any heavy rotor molecules (e.g. HCN) 
and their isotopologues. Chemical signatures of very high ionization fractions 
in dense UV-shielded gas such as low [DCO+]/[HCO+] and high [HCO+]/[CO] 
abundance ratios would be good probes of CRDRs in extreme starbursts. These 
tests, along with direct measurements of the high CO line brightness temper- 
atures expected over the areas of compact dense gas disks found in ULIRGs 
will soon be feasible as sub-arcsecond interferometric imaging capabilities and 
sensitivity at mm/submm wavelengths improve in the era of ALMA. 

Subject headings: galaxies: star formation - galaxies: starburst - galaxies: IR - 
ISM: molecules - ISM: dust - ISM: cosmic rays 
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1. Introduction 

Cosmic rays (CRs) have been established for some time as the main regulators of the 
temperature, ionization, and chemical state of dense gas cores lying deep inside the (far- 
UV)-shielded regions of molecular clouds (e.g. Goldsmith & Langer 1978; Goldsmith 2001; 
Lequeux 2004 and references therein). The association of CRs with 0,B star clusters and 
supernovae remnants (SNRs) where they are accelerated (in massive-star winds and SN 
shocks) has been recently demonstrated for the Galaxy (Binns et al. 2008), and even shown 
for individual SNRs (Acciari et al. 2009a) and star-forming (SF) regions (Abdo et al. 2010), 
while the synchrotron emission of CR electrons is a well-established marker of star-forming 
regions (e.g. Condon et al. 1990, 1991). The recent detections of 7-rays (the product of 
inelastic collisions between CR-protons and hydrogen nuclei in the ISM) from the starburst 
nuclei of M82 and NGC253 solidified the connection of CR energy density to SNRs and star 
formation activity in galaxies (Acciari et al. 2009b; Acero et al. 2009). CRs rather than far- 
UV photons have even been advocated as the dominant heating mechanism of molecular gas 
in galaxies irrespective of their star formation activity (Suchkov, Allen, & Heckman 1993 and 
references therein). These early proposals however faced two distinct problems: a) ensembles 
of standard Photon Dominated Regions (PDRs) accounted well for the global molecular 
and atomic line emission from quiescent spirals such as the Milky Way (Fixsen et al. 1999; 
Mochizuki & Nakagawa 2000) as well as starbursts (e.g. Wolfire et al. 1990; Mao et al. 2000), 
and b) a tight correlation between CO line brightness and non-thermal radio continuum 
used as evidence for CR-heating of molecular gas can be also attributed to the well-known 
far-IR/radio correlation. The latter is established by star formation powering both the far- 
IR dust continuum and the non-thermal radio emission in galaxies, with far-UV photons 
from SF sites heating the dust (and thus the gas via photoelectric heating) and reprocessed 
into the far-IR continuum (e.g. Condon & Yin 1990; Condon 1992). The tight far-IR/radio 
correlation and the close association of 0,B stars and SNRs (the sites of both far-UV photons 
and CR acceleration) down to individual CO-luminous Giant Molecular Clouds (GMCs) can 
then easily account for the observed (CO intensity) / (radio continuum) correlation without 
CR-heating of CO-bright clouds as an underlying cause (i.e. the far-IR/radio becomes a 
CO/radio correlation as CO luminosity is a good proxy for far-IR luminosity). 

It was the recent capability for sensitive observations of high-excitation high-J transi- 
tions of CO and 13 CO that "broke" the aforementioned degeneracies and definitively demon- 
strated the presence of CR-heated molecular gas in the starburst nucleus of the otherwise 
quiescent spiral galaxy NGC253 (Bradford et al. 2003; Hailey-Dunsheath et al. 2008), while 
strong evidence suggests this is also the case for the dense molecular cloud in Sgr B2 at the 
Galactic Center (Yusef-Zadeh et al. 2007). Nevertheless these amount to only ~0.1-1% 
of the total molecular gas found in typical spirals and thus do not change the paradigm of 
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far-UV photons as the main heating agent of molecular gas in IR-luminous galaxies. It must 
be noted that irrespective whether CRs heat most of the molecular gas in galaxies or not, 
they remain the ultimate regulator of its temperature and ionization fraction minima, both 
reached in its densest, UV-shielded, phase deep inside molecular clouds. 



2. The CR energy density in ULIRGs: a CR-dominated ISM 

It is the star formation rate (SFR) density rather than the total SFR that determines 
the CR energy density Ucr (eVcm -3 ) in the ambient ISM of a galaxy. Indicatively, for the 
100 pc starburst nucleus of NGC253 the latter is ~2 x 10 3 times higher than the Galactic 
value (~0.5eVcm -3 ) even though their global SFRs are similar (Accero et al. 2009). In 
the central 500 pc of the nearby starbust M82 it is UcR~500xUcR,Gai (Suchkov et al. 1993; 
Acciari et al. 2009b) while its globally averaged CR energy density remains similar to that 
of the Galaxy. Such high values of Ucr can be easily attained and surprassed over the entire 
volume of the substantial reservoirs of molecular gas (~(10 9 -10 10 ) M Q ) in Ultra Luminous 
Infrared Galaxies (hereafter ULIRGs) and their high redshift counterparts. Indeed the large 
IR luminosities of these extreme starbusts (Li R (8-1000/zm)>10 12 L ) emanate from very 
small volumes with typical IR brightnesses of 

<j IR (ULIRGs) = (10 12 - 10 13 ) (1) 

and with most such systems strongly clustering around <7ir~10 13 L q kpc -2 , (Thompson, 
Quataert, & Murray 2005). This latter value could be indicating radiation-pressure reg- 
ulated maximal starbusts (Thompson, Quataert, & Murray 2005; Thompson 2009), where 
an Eddington limit from O, B star clusters sets a maximum gas accretion rate onto SF sites 
deep inside molecular clouds via photon pressure on its concomitant dust (see also Scoville 
2004 for an earlier and simple exposition). Interestingly a similar threshold value for <tir can 
be also recovered with CRs instead of photons setting the Eddington limit (Socrates, Davis, 
& Ramirez-Ruiz 2008) as CRs are much more highly coupled to the ISM than photons. 
Either way the high <ti R values typical in ULIRGs seem to be the result of extreme star- 
bursts occuring in very compact regions, while for less IR-luminous systems (Lir<10 u L ), 
which typically have more extended star formation: <Ti R = (10 10 -10 n ) L kpc -2 (Lehnert & 
Heckman 1996). For nearby ULIRGs this compactness of their molecular gas reservoirs 
(and thus of their star-forming volumes) has been demonstrated using mm and (recently) 
submm interferometric imaging of their CO line and dust continuum, that found gaseous 
disks with D~(100-300) pc (e.g. Downes et al. 1998; Sakamoto et al. 2008; Matsushita et al. 
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2009). Moreover, even when such systems initially seem as having larger dimensions, higher- 
resolution mm/submm interferometry frequently reveals them as mergers of such compact 
nuclei, each containing the molecular gas of a gas-rich progenitor with M(H 2 )~(10 9 -10 10 ) M Q 
(Evans et al. 2002; Sakamoto et al. 2008). 

Recently ctir values similar to those typical for local ULIRGs have been found also 
for a submillimeter-selected galaxy (SMG) at z~2.3 where a unique combination of high- 
resolution submm imaging and a strong magnification by gravitational lensing made possible 
the resolution of the star-forming area of this distant ULIRG at linear scales of ~100pc 
(Swinbank et al. 2010). Finally, while the level of the contribution of an Active Galactic 
Nucleus (AGN) to the tremendeous <ti R values of such compact regions remains a matter of 
debate (e.g. Downes & Eckart 2007), it can be safely assumed that they are good order- 
of-magnitude "calorimeters" of dust-obscured SFR (e.g. in the archetypal QSO/starburst 
system Mrk231: L m =2/3(starburst) + l/3(AGN), Downes & Solomon 1998). This is further 
supported by the large number of radio SNRs found recently within the gaseous disks of 
ULIRGs using VLBI imaging (Sakamoto et al. al. 2008 and references therein; Perez- Torres 
et al. 2009), while similar AGN contributions to total IR luminosities are recovered also for 
dusty starbursts at high redshifts (Pope et al. 2008; Murphy et al. 2009). 

The IR luminosity surface density and the SN rate surface density can then be related 
using the well-established IR/radio correlation 



(Yun et al. 2001 and references therein). This relates the power uL u of the non-thermal radio 
continuum to the IR luminosity (and holds over five orders of magnitude with a dispersion 
< 80%), with the star formation rate consider as the ultimate scaling factor for both IR 
and non-thermal radio luminosities (see Thompson et al. 2006 for the latest theoretical 
background). To obtain a working value for n{y) we follow Yun et al. 2001 and use their 
value for the far-IR/radio correlation at z/=1.4GHz: (uL u ) /L FIR =2.66 x 10~ 4 x 10~ q , where 
q=2.35 (as defined by Helou et al. 1985). For (L m (8-1000/xm)/L F i R (42-120/W))=1.3 (Yun 
et al. 2001), the coefficient in Equation 2 becomes /x(z/) = 1.54, and is adopted for this study 
(with an inherent uncertainty of a factor of ~2 due to a varying L FIR /L IR in IR-luminous 
galaxies). On the other hand the non-thermal radio power depends on the SN rate fsN as 




H{v) 10 



(2) 



lO^WHz" 1 



13 KGHzjPWyr- 1 ), 



(3) 



(Condon 1992). Combining Equations 2 and 3 and assuming full concomitance of IR and 
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non-thermal radio emission over the star forming regions yields the SN rate surface density 



= 2 ' 95 " M Goto*) MGHz)r1 ' (4) 

Following Suckhov et al. (1993) the CR energy densities in ULIRGs scale with respect 
to that of the Galaxy as 



UcR Ssn / Vdiff 

~ x I 

U C R,Gal £ S N,Gal \ V wind 

where SsN,Gai=3.85xlO _5 yr _1 kpc -2 (McKee & Williams 1997), and V diff is the diffusion 
velocity at which CRs escape from quiescent disks like the Milky Way while V w i n d is the 
velocity of a SF-induced wind at which CRs are advected out of the SF regions in star- 
bursts. Typically Vdifi~10kms _1 while the maximum velocity of SF-driven galactic winds 
is V w i n d ima x~(2-3)xl0 3 kms _1 (Suchkov et a. 1993; Veilleux et al. 2005). Such starburst- 
induced high wind velocities are deduced for M82 (Seaquist et al. 1985) though usually 
V w i n d~(500-1000) kms -1 . From Equations 4 and 5, after setting V win d, m ax=3000 km s -1 , 
and /i(^)~ 1.54, a=0.8 (for z/=1.4GHz), Equation 1 yields CR energy density in compact 
starbursts (CSB) of 




Ucr.csb 

U C R. 



Gal 



4 x(10 2 - 10 4 ). (6) 



For ULIRGs as the template compact starbursts: <r IR = 10 13 L Q kpc -2 , yielding U C r,csb~ 
4 x 10 3 UcR,Gai- The latter amounts to a tremendeous boost of CR energy density, similar 
to that found recently for a single molecular cloud in the nucleus of NGC 253 (Acero et al. 
2009), and is capable of turning the massive and dense molecular clouds of ULIRGs into CR- 
dominated Regions (hereafter CRDRs) in terms of the dominant heating mechanism. For 
the extreme starbursts in ULIRGs and their high redshift counterparts these high <7ir and 
correspondingly high Ucr values can involve ~(10 9 -10 10 ) M of molecular gas mass fueling 
galaxy-wide star-forming episodes. 



2.1. CRDRs in ULIRGs: a new temperature minimum for UV-shielded gas 

For the UV-shielded, and mostly subsonic, dense gas cores deep inside molecular clouds 
the heating rate is given by 
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where Ccr(s _1 )oc U cr is the CR ionization rate per H 2 molecule. The gas cooling via gas-dust 
interaction can be expressed as 



AH - 10- 25 (i^s) 2 Ti /2 (T k - T dust ) erg cm" 3 s" 1 (8) 

(see Tielens 2005 for derivation of both expressions). The most important line cooling of 
this gas phase is through rotational lines of CO and a few other molecular species whose 
rotational ladder energy levels reaches down to AE ul /k B ~(5-10) K. Following the detailed 
study of Goldsmith (2001) we parametrize the molecular line cooling as 



1/2 / ^ n /3 

(lblw er s cnr3s_1 - ( 9 ) 

The density dependence was extracted from a fit of the parameter a in Table 2 of Goldsmith 
(2001), and reproduces the values of the A linc to within < 20% for n(H 2 ) = (10 4 -10 6 ) cm" 3 , 
which spans the density range of dense cores within GMCs. For that density range /3~3. 

The resulting gas temperature for the CR-heated gas can then be estimated from the 
equation of thermal balance 

r C R = Aline + Ag-d- (10) 

Setting the dust temperature to T dust =0 K will yield a minimum T k value while also allowing 
a simple analytic solution of Equation 10. Substituting the expressions from Equations 7, 8 
and 9 into the latter yields 

T 3 il0 + 0.526nf T^ 2 = 0.25nfc_i 7 , (11) 

where T k) i =T k /(10K), n 4 =n(H 2 )/(10 4 cnT 3 ) and C-it=Ccr/(10" 17 s" 1 ). An exact solution 
of the latter then is 



A 



line 



6 x lO" 24 



n(H 2 ) 
10 4 cm" 3 



-k,10 



0.630 



n 



1/2 



C-i7 + 0.276676n 3 V /2 - 0.526n4 /2 



2/3 



(12) 
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For CcR,Gai=5xlO _17 s _1 for the Galaxy (e.g. van der Tak & van Dishoeck 2000) and 
n(H 2 )=10 4 cm -3 the latter yields T k ~9K which is typical for UV-shielded gas immersed 
in Galactic CR energy density (e.g. Goldsmith 2001), and deduced by numerous observa- 
tions in the Galaxy (e.g. Pineda & Bensch 2007; Bergin & Tafalla 2007). For Ccr,csb=(1~ 
4) x 10 3 CcR,Gai expected for the ISM of extreme starbursts, and n(H2) = (10 4 -10 6 ) cm -3 (typ- 
ical densities for dense cores in the Galaxy e.g. Bergin & Tafalla 2007) Equation 12 yields 
T k ~(80-240) K, as the minimum possible temperature for the molecular gas in compact ex- 
treme starbursts (see Figure 1). Turbulent gas heating (e.g. Pan & Padoan 2009) that may 
"seep" down the molecular cloud hierarchical structures to the dense gas cores (although as 
typically subsonic, they are expected to have neglible such heating), any sort of mechanical 
heating of the dense gas in ULIRGs (Baan, Loenen, & Spaans 2010), or warmer dust because 
of IR light "leaking" deep inside molecular clouds can only raise this temperature range. 

The gas cores deep inside the CR-heated regions of molecular clouds in the Galaxy, 
and especially the highest density ones (~(10 5 -10 6 ) cm~ 3 ) are typically dominated by near- 
thermal motions (e.g. Bergin & Tafalla 2007). However the often large turbulent linewidths 
found in the dense gas disks of nearby ULIRGs (Downes & Solomon 1998; Sakamoto et 
al. 2008, Matsushita et al. 2009) could in principle affect even dense gas core kinematics 
and thus their line cooling function. Following Goldsmith (2001) that any such macroscopic 
motions are driven mostly by self-gravity, a new effective line cooling function would be 



which rougly quantifies the effect of increased transparency (and thus cooling power) of 
molecular line photons. An optical depth dependance of roc(dV/dR)~ 1 is adopted, where 
(dV/dR) V iR=[n(H 2 )/(10 3 cm- 3 )] 1/2 kms" 1 pc 1 is the velocity gradient of macroscopic mo- 
tions for self-gravitating gas. The equation of thermal balance and its solution then become 




(13) 



T 3 10 + 0.166n 4 T^ 



= 0.0789C-17, 



(14) 



and 



0.630 [(0.02766n^ + 0.3156C-i 7 ) 1/2 - 0.1663n 4 



2/3 



(15) 



(where we set again Td us t=0K). For CcR~(l~4)xl0 3 CcR,Gai the new thermal balance equa- 
tion yields T k (min)~(55-115) K for the molecular gas in such environments while for the 
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more extreme CSBs with CcR~10 4 CcR,Gai it is T k (min)~(145-160) K. In all cases the mini- 
mum temperatures in CRDRs remain significantly higher the minimum Tk~(8-10) K which 
is attained in the dark UV-shielded cores in the Galaxy (see Figure 1) where the initial 
conditions of star formation are set (Bergin & Tafalla 2007). 

It must be noted that such high CR energy densities may also affect molecular gas chem- 
istry, and thus the abundance of coolants such as CO. Hence better temperature estimates of 
the UV-shielded dense gas cores in CRDRs can only be provided by self-consistent solutions 
of their thermal and chemical states. Moreover activation of other cooling lines such as 01 
at 63yum (AE u i/kB~228 K) when temperatures rise significantly above 100 K in dense gas 
(n(H 2 ) > 10 5 cm- 3 ) can cap the rise of T k (min)=F(CcR) in CRDRs to ~ 150 K (Thi 2010). 
One the other hand gas temperatures can be even higher than the simple estimates provided 
by Equation 10 because of residual turbulent motion dissipation (and thus heating) remains 
possible in the dense gas reservoirs of ULIRGs, as they may resemble the tidally-stirred 
dense gas in the Galactic Center where turbulent heating remains important even at high 
densities (e.g. Stark et al 1989; Rodriguez- Fernandez et al. 2001; Giisten & Philipp 2004). 



2.2. The ionization fraction of dense gas in ULIRGs 

The large cosmic ray energy densities, besides significantly raising the minimum pos- 
sible temperature of molecular gas in the compact starbursts powering ULIRGs and their 
high redshift counterparts, they will also dramatically raise the minimum ionization frac- 
tion. Following the treatment by McKee (1989), in UV-shielded environments with negli- 
gible photoionization and CRs as the sole cause of ISM ionization, the ionization fraction 
x(e)=n c /n(H) = n e /2n(H 2 ) is given by 



x(e) = 2 x 10~V 



n ch 



gd 1 2n(H, 



1/2 



1 + 



n c h 



8n(H 2 



1/2 



+ 



8n(H 2 ) 



1/2' 



(16) 



where n ch ~500 (rg d C-17) cm" 3 is a characteristic density encapsulating the effect of cosmic 
rays and ambient metallicity on the ionization balance (r gd is the normalized gas/dust ra- 
tio r gd =[(G/D)/100] with G/D(gas-to-dust mass)=100 assumed for Solar metallicities; e.g. 
Knapp & Kerr 1974; Aannestad & Purcell 1973). 

For the Galaxy where CcR,Gai=5x 10~ 17 s _1 it is n ch =2.5x 10 3 cm~ 3 , and for a typical 
dense core density of n(H 2 ) = 10 5 cur 3 Equation 16 yields x(e)~2.4x 10 -8 , consistent with 
the typical range in the Galaxy: 5 x 10~ 9 < x(e) < 1.5 x 10~ 7 (e.g. Langer 1985). For 
the much larger CR ionization rates CcR,csB=10 3 x£cR,Gai expected in the ISM of compact 
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starbursts: n ch =2.5x 10 6 cm -3 and x(e)~4xl0~ 6 . The latter is one to two orders of mag- 
nitude larger than the typical range and ~4 times higher than the highest value measured 
for dense UV-shielded cores anywhere in the Galaxy (Caselli et al. 1998). In the classic 
photoinization-regulated star formation scenario (McKee 1989) such high ionization frac- 
tions will be capable of turning even very dense gas in the UV-shielded/CR- ionized regions 
of molecular clouds subcritical (i.e. M core /M$<l where M$=0.12$/G 1 / 2 , and $ is the mag- 
netic field flux threading the molecular core, Mouschovias & Spitzer 1976), and thus halt 
their gravitational collapse until a now much slower ambiplolar diffusion allows it. 



3. Important consequences and some key observational tests 

The CR-permeated molecular gas in compact extreme starbursts is more than the mere 
sum of individual star-forming regions and their localized dense PDRs. The latter would 
leave most of the dense gas settle to a cold state since for the high gas densities found in 
ULIRGs, far-UV field intensities will be reduced by factors of ~10 4 over distances of <0.1 pc. 
On the other hand, by dramatically altering the thermal and ionization state of dense UV- 
shielded gas cores inside molecular clouds the large CR energy densities in extreme starbusts 
significantly alter the initial conditions for star formation in such systems. Indeed it is in the 
UV-shielded dense gas cores where these initial conditions are set, and the large temperatures 
expected for these cores throughout CRDRs invalidates the main arguments about an almost 
constant characteristic mass of young stars in most ISM environments including starbursts 
(Elmegreen, Klessen, & Wilson 2008). The latter study did not consider the effects of CRs, 
and as a result found that UV-shielded dense gas remains cold (T k ~10K) even in starbursts. 

The much larger ionization fractions that can now be reached deep inside molecular 
clouds can in principle: a) keep the magnetic field lines strongly "threaded" onto molecular 
gas at much higher densities, and b) as a result render much of its mass incapable of star- 
formation, at least in the simple photoinization-regulated SF scenario. These effects stem 
from the now much longer ambipolar diffusion timescale 



r AD = 3.2 x 10 7 r gd 



n ch 



2n(H 2 ) 



1/2 



1 + 



n c h 



8n(H 2 ) 



1/2 



+ 



rich 



(~8n(H 2 ) 



1/2' 



yrs, 



;i7) 



needed for a CR-ionized dense core with density n(H 2 ) to lose magnetic flux and collapse 
(McKee 1989). For the CR ionization rates of (cR=fewx (10 3 -10 4 )(cR,Gai expected in CRDRs 
it would be n c h=2.5x 10 6 ~ 7 cm -3 , thus for a typical dense core density of n(H 2 ) = 10 5 cm~ 3 , 
t ad~3x (10 8 -10 10 ) yrs (see Figure 2). In photoinization-regulated star formation this is a 
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lower limit on the gas consumption timescale by the latter process, and in CRDRs it is 
clearly already much longer (by up to two orders of magnitude) than the typical consump- 
tion timescale of molecular gas reservoirs of LIRGs (Figure 2). For the more vigorously 
star-forming ULIRGs, and considering only their dense HCN-bright molecular gas phase 
as the true SF fuel (e.g. Gao & Solomon 2004): r cons =M(n>10 4 cm- 3 )/SFR~10 7 yrs, and 
TAD/T C o„s(n>10 4 cm- 3 )~10 2 -10 3 . 

This certainly disfavors a simple quasi-static photoionization-regulated SF scenario of 
B-field lines slowly slipping from stationary dense gas cores which then proceed to star 
formation. Such a simple picture is expected to be modified anyway by the presence of 
MHD turbulence which can accelarate the ambipolar diffusion process (Heitsch et al. 2004), 
especially in the very turbulent molecular gas of ULIRGs. For such systems the magnetic 
fields can be strong, in possible equipartition with highly turbulent gas (Thompson et al. 
2006), and thus dynamically important and strongly co-evolving with the molecular gas. 
It may well be that the ability of CRDRs to maintain high ionization fractions in very 
dense UV-shielded molecular gas, and thus retain a strong coupling of the magnetic fields 
onto the bulk of its mass in ULIRGs (where (n(H 2 ))>10 4 cm -3 ), is what allows a quick 
establishment of equipartition between magnetic fields and turbulent gas motions in their 
self-gravitating disks. MHD simulations for gas clouds immersed in the very intense CR 
energy density backrounds expected in CRDRs will be key in addressing these issues, and 
investigate whether turbulence can still accelerate ambipolar diffusion in molecular clouds 
of such high ionization fractions and lead back to TAD<Tcons f° r the dense gas in ULIRGs. 



3.1. A new set of initial conditions for star formation in ULIRGs 

In the high-extinction ISM of compact extreme starbursts, CRDRs make the dramatic 
change of star-formation initial conditions an imperative of high- density star formation. The 
latter occurs irrespective whether the bulk of the molecular gas in CRDRs of ULIRGs is 
heated by CRs or not as the latter will now raise the temperature minimum of pre-stellar 
UV-shielded dense gas cores throught their star-forming volumes by a large factor. The 
consequences of this new imperative for star-formation in extreme starburst environments 
remain invariant irrespective whether gravitational instability in turbulent molecular gas 
(e.g Klessen 2004; Jappsen et al. 2005; Bonnell, Clarke, & Bate 2006) or ambipolar diffusion 
of magnetic field lines from dissipated dense cores followed by their gravitational collapse 
(e.g. Mouschovias & Spitzer 1976, McKee 1989) drive star formation in galaxies. In both 
schemes the dense UV-shielded cores of molecular clouds is where the initial conditions of 
star formation are trully set (see Elmegreen 2007; Ballesteros-Paredes & Hartmann 2007 for 
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recent excellent reviews). It must also be noted that this gas phase is different from the 
PDR-dominated gas around star-forming sites that dominates the global molecular line and 
dust continuum spectral energy distributions observed for ULIRGs and often (erroneously) 
used to set the star formation initial conditions in starbursts (e.g. Klessen et al. 2007). 

The effects on the characteristic mass scale of young stars M^. h ; and thus on the stellar 
IMF (Elmegreen et al. 2008) for the ISM in CRDRs are explored in detail in a forthcoming 
paper (Papadopoulos et al. 2010). It nevertheless worths pointing out that the large boost 
of Tk(min) in CRDRs in an (almost) extinction-free manner across a large range of densities 
in molecular clouds cannot but have fundamental consequences on M^ h and the emergent 
stellar IMF. Indicatively for a Tk(min) in UV-shielded cores boosted by a factor of 10, the 

Jeans mass MjocT^ 2 n(H 2 ) -1 / 2 rises by a factor of ~32! (over an identical density range), 

(*) 

and almost certainly raises M^; and the characteristic mass scale of the stellar IMF. Inter- 
estingly similar effects can occur also in AGN-induced X-ray Dominated Regions (XDRs; 
Schleicher, Spaans, & Klessen 2010) since X-rays just as CRs (and unlike far-UV photons) 
can volumetrically heat large columns of molecular gas while experiencing very little extinc- 
tion. The effect of X-rays on the Jeans mass of dense cores and the IMF has been recently 
proposed for a powerful distant QSO (Bradford et al. 2009) and it may represent a neglected 
but important AGN feedback factor on its circumnuclear star formation. 



3.2. CRDRs: observational tests 

The molecular line diagnostics of starburst-induced CRDRs and AGN-originating XDRs 
can be to a large degree degenerate when galaxies host both power sources. This has been 
noticed in earlier comparative studies of XDRs and regions with higher Ucr values (though 
only up to 100xUcR,Gai) found that only carefully chosen line ratios can distinguish between 
them (Meijerink, Spaans, & Israel 2006). The still larger CR energy densities expected in 
CRDRs of ULIRGs will further compound these difficulties. 

Provided that a powerful X-ray luminous AGN heating up the bulk of the molecular 
gas in its host galaxy can be somehow excluded (e.g. via hard X-ray observations), any set 
of molecular lines and ratios that can strongly constrain the temperature of the dense gas 
UV-shielded phase (n(H 2 )>10 4 cm" 3 ) in ULIRGs will be valuable. Indeed, given that in the 
hierachical structures of typical molecular clouds the dense gas regions: a) lie well inside 
much larger ones that strongly attenuate far-UV light and b) cool strongly via molecular 
line emission because of their high densities, then any evidence for high temperatures for the 
dense gas phase would be a indicator of strong CR-heating. In that regard observations of 
high-J CO lines such as J=6-5 and its 13 CO isotopologue have already been proven excellent 
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in revealing CR- heated rather than UV/photoelectrically- heated molecular gas in galactic 
nuclei (Hailey-Dunsheath et al. 2008). Irrespective of the particular set of rotational lines 
used as a "thermometer" of the dense gas, three general requirements must be met: a) all lines 
must have high critical densities (n crit >10 4 cm" 3 ), b) have widely separated Ej +1 j/1cb factors, 
and c) the J-corresponding lines of at least one rare isotopologue must also be observed (e.g. 
12 CO and 13 CO, or C 32 S and C 34 S, etc). The first two ensure probing of the dense gas phase 
while maintaining good T k -sensitivity, and the last one is necessary for reducing well-known 
degeneracies when modeling only transitions of the most abundant isotopologue which often 
have significant optical depths. Molecular lines with n crit >10 5 cm -3 (e.g. HCN, CS rotational 
transitions) are particularly valuable since, aside from emanating from gas well within typical 
pre-stellar cores, they trace a phase whose kinematic state is either dictated by self-gravity 
(e.g Goldsmith 2001), or has fully dissipated to thermal motions. This constrains the line 
formation mechanism (and can be used to remove degeneracies of radiative transfer models 
e.g. see Greve et al. 2009), but most importantly it reduces the possibility of residual 
mechanical heating of the dense gas (Loenen et al. 2008) that could mask as CR-heating (as 
both can heat gas volumetrically) . 

A brief example of such diagnostics can be provided using a Large Velocity Gradient 
(LVG) code (e.g. Richardson 1985) to compute relative line intensities for dense gas with 
T k =(10-15)K, and n(H 2 ) = 10 5 cm" 3 , and T k =(100-150) K at the same density. In both 
cases a gas velocity gradient due to self-gravity is assumed. For the cold gas, the CO 
and 13 CO (J+l-J)/(3-2) brightness temperature ratios are: R(j+i,j)/32 ^0.7 for J+l>4 
(J=l-0, and 2-1 are not considered because they can have substantial contributions from 
a diffuse non self-gravitating phase). For the warm gas R(j+i 5 j)/32~0.9-0.95 for J+ 1=4-6, 
while the corresponding ratios for 13 CO are ~1 1.3. Similar diagnostics but using multi-J line 
emission from rarer molecules (and their isotopologues) with much larger dipole moments 
(n cr i t >10 5 cm" 3 ) such as HCN and H 13 CN can even better constrain the temperature of 
dense gas in galaxies. However their much fainter emission (e.g. HCN lines are ~5-30 times 
fainter than those of CO) allows their use only for the brightest nearby starburst nuclei (e.g. 
see Jackson et al. 1995 for an early pioneering effort), and only ALMA will enable such 
diagnostics for a large number of galaxies in the local and distant Universe. 

Strong thermo-chemical effects induced by large U C r values and their ability to volu- 
metrically warm large amounts of dense molecular gas can also provide valuable diagnostic of 
the temperatures deep inside dense gas cores. For example global HNC/HCN J=l-0 bright- 
ness temperature ratios of R^?, HCN =0.5-1.0 in galaxies can be attributed to an ensemble 

of PDRs but ratios Rhnc/hcn^-5 cannot, and imply T k >100K (Loenen et al. 2008). Such 
low ratios are indeed found in ULIRGs but are attributed to turbulent gas heating by SNRs 
in dense molecular environments (Loenen 2009) that substantially warms the gas in the 
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absense of photons to T k >100K necessary for the HNC+H— *-HCN+H reaction to proceed 
efficiently and convert HNC to HCN (e.g. Schilke et al. f992). Unlike the Galaxy where 
turbulent gas heating has subsided in the dense subsonic gas cores that would emitt in these 
transitions this may not be so in the ISM of ULlRGs, making CR and mechanical heating 
difficult to distinguish. As mentioned earlier, moving the CRDR molecular line diagnostic 
to tracers with ever increasing critical densities (and definitely with n cr ; t >10 5 cm~ 3 ) may 
be the only way of "breaking" this degeneracy as turbulent heating is expected to become 
progressively weaker in higher density cores (after all star formation is expected to proceed 
in dense gas cores whose turbulence has fully dissipated, Larson 2005). Finally, chemical 
signatures that can uniquely trace the high Ucr values in CRDRs using the high ionization 
fractions expected for their dense gas cores are particularly valuable. A very sensitive such 
probe of Ccr and x(e) in dense UV-shielded gas is provided by Rd=[DCO + ]/[HCO + ] versus 
Rh=[HCO + ]/[CO] abundance ratio diagrams where the high CR energy densities and result- 
ing ionization fractions in CRDRs would correspond to very low Rd and very high Rh values 
(Caselli et al. 1998). A dedicated study of these issues that explores unique diagnostic of 
CRDRs is now in progress (Meijerink et al. 2010). 

Most of the aforementioned tests will acquire additional diagnostic power once ALMA 
will be able to conduct them at high, sub-arcsecond resolution in nearby ULIRGs. It will then 
be possible to discern whether very warm dense gas with high ionization fractions is localized 
around a point-like source, the AGN (and is thus due to XDRs), or is well-distributed over 
the star-forming area of compact systems and is thus starburst-related (see Schleicher et 
al. 2010 for a recent such study). High resolution imaging can also directly measure the 
high brightness temperatures expected for all optically thick and thermalized molecular 
lines (where Tex De ^~Tk) in CRDRs. Radiative transfer models of emergent CO line emission 
from a self-gravitating gas phase with n(H 2 )=10 5 cm~ 3 and Tk=(100-150) K typically yield 
T\,(CO)~(80-140) K for J=l-0 up to J=7-6, and imaging at a linear resolution of AL<100pc 
(the diameter of gaseous disks in ULIRGs) would be adequate to directly measure them (see 
Sakamoto et al. 2008; and Matshushita 2009 for early examples). For z<0.05 (which would 
include a large number of ULIRGs) such linear resolutions correspond to angular resolutions 
of 6'b~0.1 // (for a flat A-dominated cosmology with H = 71kms~ 1 Mpc~ 1 and f2 m = 0.27) 
which would be certainly possible with ALMA. 
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4. Conclusions 

The main conclusions of this work are three, namely 

• In the high star-formation density environments of ULIRGs cosmic ray energy densities 
Ucr will be enhanced by a tremendeous factor of ~fewx (10 3 -10 4 )UcR,Gai- These will 
permeate the large molecular gas reservoirs of such systems, likely turning them into 
CR-dominated regions (CRDRs) where CRs not far-UV photons regulate the thermal 
and ionization state of the bulk of their typicaly very dense molecular gas. 

• Irrespective whether the tremendeous CR energy densities in the compact starburst 
regions of ULIRGs provide the dominant heating for most of their molecular gas or 
not, they will dramatically raise the gas temperature and ionization fraction minimum 
values possible in their ISM, which are typically attained in UV-shielded dense gas 
cores where the star-formation initial conditions are set. 

• The new and very different initial conditions for star formation in CRDRs are an 
imperative for all high-density star-formation events, and will almost certainly raise 
the characteristic mass of young stars (and thus the stellar IMF mass scale) during 
such events throughout the Universe. 

Sensitive observations of key molecular lines with high critical densities (n crit >10 5 cm -3 ) 
as well as high resolution mm/submm interferometric imaging will be key tools in uncovering 
the high temperatures of the dense UV-shielded gas in CRDRs expected in compact extreme 
starbursts. Tracers of the very high ionization fractions expected for their dense gas can 
provide an independent assesment of their presense. All these tests will become possible 
with ALMA for large numbers of ULIRGs in the local Universe. 

It is a pleasure to acknowledge many informative discussions with Drs Wing-Fai Thi 
and Francesco Miniati. Discussions with Dr Marco Spaans regarding XDR/CDR potential 
diagnostics are also greatfuly acknowledged. Finally I would like to thank the referee for 
his/her extensive set of comments that helped helped clarify several issues and significantly 
widen the scope of this work. 
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Fig. 1. — The temperatures of UV-shielded gas in CR-dominated Regions (CRDRs) of 
compact extreme starbursts in ULIRGs (marked by shaded area) for cores with densities 
n(H 2 ) = (10 4 ,10 5 ,10 6 ) cm -3 and near thermal motions (left), or dictated by gas self-gravity 
(right) that yield larger than thermal linewidths and stronger line cooling (see section 2.1). 
The horizontal line marks the typical temperature of this gas phase in the Galaxy. The 
arrow marks the adopted Galactic CR ionization rate of CcR,Gai=5 x 10~ 17 s _1 . 
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Fig. 2. — The ambipolar diffusion timescale of UV-shielded gas in CR-dominated Regions 
(CRDRs) of compact extreme starbursts in ULIRGs (marked by shaded area) for cores 
with densities n(H 2 ) = (10 4 ,10 5 ,10 6 ) cm -3 estimated from Equation 17. The lower horizon- 
tal line marked by TAD,Gai is the ambipolar diffusion timescale for UV-shielded cores with 
n(H 2 ) = 10 5 cm -3 in the Galaxy for CcR,Gai=5 x 10~ 17 s _1 , and the higher line marks the gas 
consumption timescale of the entire molecular gas reservoir of a typical LIRG (L IR ~10 11 L ). 
For the more vigorously star-forming ULIRGs and their HCN-bright dense gas phase this 
gas consumption timescale can be an order of magnitude less (section 3). 



